Effects of Multi-Surface Modification on Curie temperature of 

ferroelectric films 



Xiao-Guang Wang* 

Laboratory of optical physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100080, P. 
R.China 

Shao-Hua Pan and Guo-Zhen Yang 
China Center of Advanced Science and Technology (World Laboratory), P.O.Box 8730, Beijing 
100080,P. R. China and Institute of Physics, Chinese Academy of Sciences, Beijing 100080, P.R. 
China 



Abstract. Within the framework of mean field theory, we study the effects of multi-surface 
modification on Curie temperature of ferroelectric films using the transverse Ising model. The general 
nonlinear equations for Curie temperature of multi-surface ferroelectric films with arbitrary exchange 
Q\ constants and transverse fields are derived by the transfer matrix method. As an example, we consider 

a (l,n,l) film consisting of / top surface layers, n bulk layers and / bottom surface layers. Two 

types of surface modifications, modifications of a surface exchange constant and a surface transverse 
^ field are taken into account. The dependence of Curie temperature on the surface layer numbers, bulk 

Q layer numbers, surface exchange constants , surface transverse fields and bulk transverse fields is 

discussed. 

1 Introduction 

© Size effects are expected to be particular important in materials that undergo phase transitions, 

such as ferromagnets, superconductors and ferroelectrics[l-3]. In recent years, there has been 
considerable interest in the effects of size and surface on the ferroelectric phase transition due to the 
^ developments of ferroelectric films and composite materials [4-20]. Nonvolatile memories and 

— , other applications have made ferroelectric films a subject of great practical importance [21]. 

An early study on KDP fine particles embedded in an insulating medium showed no ferroelectric 
£h phase transition if their size is less than 150 nm [22], while Anliker et al. demonstrated that the 

I transition temperature of BaTiOs increases as the particle size decreases [23]. Controversial results 

obtained from BaTi0 3 [24] fine particles and PbTi0 3 fine particles[25] show that the transition 
g temperature decrease with decrease in grain size. 

q Theoretically, much work has been done on the Laudau phenomenological theory of phase 

transition in ferroelectric films [6,8,10,1 1]. On the microscopic level, the pseudo-spin theory based 
on the transverse Ising model (TIM) was first introduced by de Gennes [26] to describe the phase 
transition of hydrogen-bonded ferroelectrics. The model has been applied to many other systems 
such as ferromagnetic systems [27] and Jahn-teller systems [28]. Surface exchange constant and 
surface layer thickness dependence of Curie temperature of ultra-thin ferroelectric films is studied 
numerically by considering two types of surface modification of the exchange constants [29]. Curie 
temperature for ferroelectric films of arbitrary thickness has been considered and modification of 
both surface exchange constants and surface transverse field is included [30,31]. The effects of few 
surface layers ( / = 1,2 ) of ferroelectric (/,«,/) films on Curie temperature are studied in previous 

works. Little theoretical work has been devoted to the effects of multi-surface modification ( / > 3 ) 
on the Curie temperature. In practice, there exist multiple surface layers in ferroelectric films 
fabricated by present experimental techniques. 
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In this article, we study the ferroelectric (l,n,l) films of arbitrary surface layers / and bulk 

layers n and take into account two types of surface modification, modification of surface exchange 

constant / and surface transverse field Q. s . By using the transfer matrix method introduced for 

ferromagnetic superlattices[32], explicit and general nonlinear equation for Curie temperature is 
obtained. We give the equation in the following section. The dependence of Curie temperature on 
surface layer numbers, bulk layer numbers, surface exchange constants and surface transverse fields 
are discussed in section 3 and the conclusion is given in the last section. 



2. Trans verse Ising model and transfer matrix method 

We start with the TIM[26,29,30] 

h = --Y Y j s s -Y a s: , 



(1) 



where S*,S* are the X and z components of the pseudo-spin, (i, j) are plane indices and (r, r') 
are different sites of the planes, J ~ denote the exchange constants. We assume that the transverse 
field Q. t is dependent only on layer index and consider the interaction between neighboring sites. 
The spin average (S t ^ , obtained from the mean field theory [29,30] 



^— tanh(l H, \/2k B T) , 



2\H;\ 



(2) 



where H ^,0,^ y^*)) is the mean field acting on the ith spin, , k B is the Boltzman 



constant and T is the temperature. 

At a temperature close and below the Curie temperature, (S* ^ and (S^ ^ are small, I H ( l== 
equation (2) can be approximated as 



(S') = (l/2)taiih(Q I ./2Jfc Jl r) (3) 
(S; ) = [z J u (S; ) + Z (/ U+1 (S' M ) + /,, M (SU )) J[l/ 2n, )tanh(n,. / 2k B T) (4) 

Here Z and z are the numbers of nearest neighbors in a certain plane and between successive planes 
respectively. By defining 

K i =z G J ii lzJ,K i . + ,=J. M IJ,K i ._ x =J i ._JJ,m l =(SI) 

% i = (2Q, /z7)coth(Q,. /2k B T) , (5) 

equation(4) can be written as 

(r, - K l )m, - K i M m M - K^m^ = . (6) 

Here / is the interlayer exchange constant in the bulk. 

Let us rewrite the above equation in matrix form in analogy with Ref.[32] 



V ' J 



= AT 



m. 



with M i as the transfer matrix defined by 



(7) 



2 



AT.. 



1 

V 



/ 



(8) 



We assume that the ferroelectric film contains N + 1 layers with layer indices i = 0,1,2. ..N . 
From equation (7), we get 



= R 



m. 



(9) 



where R = M N1 M 2 M l represents successive multiplication of the transfer matrices. 

For an ideal film system, there exists symmetry in the direction perpendicular to the surface, 
which allows us to write m ; = m N _ i . Then, the following nonlinear equations for Curie temperature 
can be obtained from equations (6) and (9) as 

tf n [(T -* )/*o,if +(*i2 -R 21 )y o -K )/K ,]-R 22 =0. (10) 
The above equation is the general equation for Curie temperature for arbitrary exchange constants 
7„ and transverse field Q. . 



For a uniform system of iV + 1 layers with J tJ = J and Q. i = Q., equation (10) reduces to 

R x , (t - z Q I z) 2 + (R u - R 2l )(t - z 1 z) - R 22 = , 

where I = (2Q. I zJ ) COth(^ / 2k B T) . The total transfer matrix 



R = 



'T-Z /Z -O 



1 

The Nth power of a unimodular 2x2 matrix A can be simplified as [28] 



J 



A N =U N A-U N J, 

where / is the 2x2 unit matrix and U N = — — X_ ) . A + and X_ are the two 

eigenvalues of the matrix A . 

From equations (12) and (13), we obtain the matrix elements of R as follows 

=(,1- z l z)U N _ 1 - U N _ z , 
R\2 —~U N _ l ,R 2l =C/ Ar _ 1 , 

^22 — _ ^ N-2 ■ 



(ID 



(12) 



(13) 



(14) 



Here, X ± = [T — Z Q I Z ± sJ(T— Z / z) 2 —4 ] / 2 . From the definition of U N and relations 
A + + X_ = l — Zq I Z and A + X_ = 1 , we have the recursion relations 

u N =a + +x_)u N _,-u N _ 2 , 
u N+2 =a + +x_) 2 u N -2a + +x_)u N _ l+ u N _ 2 

Substituting equation (14) to equation (11) and using the above relations, we reduce the equation (11) 
to its simplest form 

U N+2 =0 (17) 

which is identical with the result of reference [26]. 



(15) 
(16) 



U N+2 can be written as 



U N+2 =sin[(N + !)</>]/ sin</) 



(18) 



for (T - Z Q I Z) 2 < 4. Here (/) = arCCOS[(T - Z Q /z)/2]. For (T - Z I zf > 4, 
(j) becomes i 6 , and the trigonometric functions become hyperbolic functions of 6 . 
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Equation (17) gives 

(2Q/ Z J )coth(Q /2k B T) = z /z + 2cos[7i/(N + 2)]. (19) 

In the limit — > 00 the bulk Curie temperature T g is determined by 

2£lco\h(£l/2k B T B ) = n Q J , (20) 

where n Q = Zq + 2z is the total number of neighbours. It can be seen that equation (20) has a 
solution only if n Q J >2£X 



3. Curie temperature of ferroelectric films with multi-surface modification 

As an example, we consider a (/,«,/) film consisting of / top surface layers, n bulk layers and 

/ bottom surface layers and study three models. In model I, we assume that pseudo-spins lie on a 

simple cubic lattice and that the transverse field Q. i = Q. . The coupling strength in a surface layer is 

denoted by / s and that in a bulk layers or between successive layers is denoted by / . Then the 
transfer matrix (equation(8)) reduces to two different types of matrices 









(y -r 


p = 




Q = 








v 1 °j 



(21) 



and the total transfer matrix 



R = pllQnpl-l (22) 

Here X = t — 47 s I J and Y — X — 4 . The nonlinear equation (10) for Curie temperature thus 
reduces to 

R U X 2 +(R l2 -R 21 )X -R 22 = 0. (23) 
The matrix elements of R can be obtained explicitly from equations(13) as 
*n = V l l l W n X 2 Y-V l 2 _ l W n _ l X 2 -IV^V^XY + 

2(V l _ 1 V l _ 2 W n _ 1 -V 2 _ x W n )X +V 2 2 WJ +V 2 1 W n _ 1 + 2V l _ x V l _ 2 W n -V 2 _ 2 W n _ x , 
R u =-R 2l =-V l 2 _ 1 W n XY + (V l \W n _ 1 +V l _ l V l _ 2 WJX +V l _ l V l _ 2 W n Y 

(24) 

+ V 2 l W n -2V l _ l V l _ 2 W n _ l -V 2 2 W n 
R 22 = -v^wj + V^W^ + 2V l _ l V l _ 2 W„ - v^w^ . 

Here V t = (< - (X>_ ) l{(X + - Ot_ ), W„ = (J3" + - J3"_ ) /(fi + -J3_), 

a ± =(X±^X 2 -4)/2 and J3 ± =(F±VF 2 -4)/2. 

For the special case of / = 1 , equation (20) reduces to 

W n+l X 2 -2W n X +W n _, =0 (25) 

which is the same as that of Ref.[30] for the structure with one top and one bottom surface layer. In 
obtaining equation (25), we have used equation (15). 

Equation (23) can be numerically solved for arbitrary / and n .Figure 1 shows the reduced 
Curie temperature t c = k B T c I J as a function of / / J for films (l,n,l) with fixed bulk layer 
numbers n = 1 but different surface layer numbers / . Curie temperature increases with increasing / 
for larger J ' I J , while it is almost independent of / when J s I J is sufficiently small. Figure 2 
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shows t c vs. J s I J for films (3, n,3) with different bulk layer numbers.The four curves n merge 
into one curve for larger J s I J , i.e. , the Curie temperature are insensitive to the bulk layer numbers 
n in this case. On the contrary, Curie temperature depends strongly on number n for 
smaller J s I J (say J s I J <1) and t c increases as n increases. 

Only one type of surface modification, the modification of surface exchange constant J s is taken into 
account in the above discussion. Next we study model II including the modification of surface 
transverse field. Let Q. s and Q. denote transverse field in surface and bulk layers respectively. The 
nonlinear equations for Curie temperature is still the equations (23) and (24) except for the change of 
the element X from X - 4 J s I J to X s - 4J s I J ,where X s = (2Q s / 7)coth(Q s / 2k B T) . In 

figure 3, we give the Curie temperature via Q. s I J . The Curie temperature decreases as 

Q. s I J increases. The transverse field causes a reduction in the Curie temperatures of the film. 

Now we consider a more realistic model, model III, in which we assume /•• = J . for both sites 

u ■> 

in surface layers, J ~ — J for both sites in bulk layers, and J ^ — J BS = J ~J between 

interface layers. The transverse fields for sites in the surface layers, for sites 

in the bulk layers, and £1 { — Q* BS — -\j£l s £l in the interface layers. The total transfer matrix 
model III becomes 



where 



M s = 



tj/j - i -i) 



M B = 



1 

t- ^ -i) 
1 o 



o 



v 



R = M<- 2 M SB M'- 2 M BS M l s -\ 



BS 4J S /J)J/J BS ~J JS /J -(t bs -4)J s /J 
(T BS -4J s /J)J/J BS 



BS 



J S 1 ' ^BS 



(26) 



M SB =M lMt M x ^ x = 



f (?BS -WbS IJ)J 2 Wbs) -Jbs'Js ~(?BS -V, 'J)J 2 HUBS) 



(t bs -4)J/J 



BS 



-J/J 



BS 



(27) 

and T BS = 2Q. BS I J COth(^ B5 / 2K g T) . For film (2,2,2) , the transfer matrix 

reduces to M SB M BS . From equations (10) and (27), we can numerically calculate the Curie 

temperature of (l,tl,l) films. In order to compare with the results of model II, we consider the 

Curie temperature vs. £1 s / / and choose parameters the same as those in figure 3. The results are 

given in figure 4. For small £1 s / / (say £1 s / / = 0. 1 ), the Curie temperatures in model III are 
larger than those in model II. It is clear that the results are qualitatively consistent with those of model 
II, while they are qualitatively different. For large £l s I J (say £l s I J = 2.5 ) and 

J s I J = 0. 1 , the Curie temperatures of model III are smaller than those in model II. It can be seen 

from both the figure 3 and figure 4 that the Curie temperature increases as J I J increases. Figure 



5 gives the Curie temperature vs. surface layer number / . The solid line in the figure corresponds to 
the Curie temperature tji of bulk system with J ' ~ — J J and Q.^ — Q. s . Let 

£1 — > £1 s , / — > / s in equation(20), the temperature can be calculated. As seen from the 
figure, the Curie temperature increases with the increase of surface layer number and approach 
asymptotically to the bulk temperature t c . 

4. Conclusions 

By making use of the transfer matrix method, an analytical and general nonlinear equation for 
Curie temperature of ferroelectric films with arbitrary exchange constant J~ and transverse field D. i 

has been derived based on the TIM. The dependence of Curie temperature on surface layer numbers, 
bulk layer numbers, surface exchange constants, surface transverse fields and bulk transverse fields 
can be easily studied by the nonlinear equation. As an example, we study three models and investigate 
the effects of multi-surface modification on Curie temperature of symmetric ferroelectric films 
(l,n,l) . In addition, the method proposed here can be applied to phase transition of both infinite and 

finite ferroelectric superlattices. It has a favorable point because of the simplicity of the numerical 
calculations. 
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Figure Captions: 

Figure 1 The reduced Curie temperature t c vs. / / J in model I for films (1,1,1) with 1 = 1 (solid 

line); / = 2 (dotted line), / = 3 (dashed line), and / = 10 (dot-dashed line). The parameter 
Q/J = 2.0. 

Figure 2 The same as figure 1, but for films (3,n,3) with n — 2 (solid line), n — 5 (dotted line), 
n = 10 (dashed line) and n = 20 (dot-dashed line). The parameter Q.I J = 2.5 . 
Figure 3 The reduced Curie temperature vs. Q. s I J in model II for films (3,2,3) with 

J s I J =0.1(solid line), J s I J =0.5(dotted line), J s I J =1.0(dashed line) and J s I J =1.5(dot- 

dashed line). The parameter Q.I J = 2.5 . 
Figure 4 The same as figure 3, but for model III. 

Figure 5 The reduced Curie temperature vs. / in model III. The parameters 
a s =2.0,a = 2.5,J s =1.2 and 71=5. 
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Figure 1 
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